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1.1 Introduction 

A main obstacle for the experimental realization of a quantum computer is the unavoidable 
coupling of the qubits to external degrees of freedom and the decoherence caused in that 
way. A possible solution of this problem are error correcting codes. These, however, require 
redundant coding and, thus, a considerably higher algorithmic effort. 

Yet another route to minimize decoherence is provided by the use of time-dependent con- 
trol fields. Such external fields influence the coherent and the dissipative behavior of a quan- 
tum system and can extend coherence times significantly. One example is the stabilization of 
a coherent superposition in a bistable potential by coupling the system to an external dipole 
field The fact that a driving field reduces the effective level splitting and therefore de- 

celerates the coherent dynamics as well as the dissipative time evolution is here of cruical 
influence. A qubit is usually represented by two distinguished levels of a more complex quan- 
tum system and, thus, a driving field may also excite the system to levels outside the doublet 
that forms the qubit, i.e., cause so-called leakage. While a small leakage itself may be tol- 
erable for the coherent dynamics, its influence on the quantum coherence of the system may 
be even more drastic. We demonstrate in this article that in a drivien qubit resonances with 
higher states, which are often ignored, may in fact enhance decoherence substantially. 

A related phemomenon has been found in the context of dissipative chaotic tunneling near 
singlet-doublet crossings where the influence of so-called chaotic levels yields an enhanced 
loss of coherence [p 



1.2 The model and its symmetries 

We consider as a working model the quartic double well with a spatially homogeneous driving 
force, harmonic in time. It is defined by the Hamiltonian 

2 1 2 4 

Hit) = -,mujlx^ + "^^^x^ + Sx cos(m). (1.1) 

^ ' Ira 4 " 64£;b 

The potential term of the static bistable Hamiltonian, possesses two minima at a; = 

ibxo, xq = {^E^jmuj^Y^'^ , separated by a barrier of height (cf. Fig. The parameter 
loq denotes the (angular) frequency of small oscillations near the bottom of each well. Thus, 
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Figure 1.1: Sk etch of the double well poten- 
tial in Eq. ([li]) for D = E-b/Tujo = 2. The 
horizontal lines mark the eigenenergies in the 
absence of the driving; the levels below the bar- 
rier come in doublets. 



the energy spectrum consists of approximately D = E-q j hujQ doublets below the baiTier and 
singlets which lie above. As a dimensionless measure for the driving strength we use F = 



The Hamiltonian (1.1 ) is T-periodic, with T = 2tt/D,. As a consequence of this discrete 
time-translational invariance of H{x, p; t), the relevant generator of the quantum dynamics is 
the one-period propagator [T 



(1.2) 



C/(T,0) =Texp l^^i^ dtHuwit)^, 



where T denotes time ordering. According to the Floquet theorem, the Floquet states of the 
system are the eigenstates of U{T, 0). They can be written in the form 



(1.3) 



with 



\Mt + T)) = \Mt))- 

Expanded in these Floquet states, the propagator of the driven system reads 

C/(t,i') = Ee""°^*"*''^''l'^"(i))(0c.(i')l- 



(1.4) 



The associated eigenphases Cq,, referred to as quasienergies, come in classes, ea.k = ea + khil, 
k = Q, ±1, ±2, .... This is suggested by a Fourier expansion of the \(j)a{t)). 
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(1.5) 



1.3 Coherent tunneling 



3 



The index k counts the number of quanta in the driving field. Otherwise, the members of a 
class a are physically equivalent. Therefore, the quasienergy spectrum can be reduced to a 
single "Brillouin zone", —Ml/2 < e < hfl/2. 

Since the quasienergies have the character of phases, they can be ordered only locally, not 
globally. A quantity that is defined on the full real axis and therefore does allow for a complete 
ordering, is the mean energy [0,p|] 



TJo 



E^ = ^ I dt i/DwW IV-a W> (1-6) 



It is related to the corresponding quasienergy by 



T 

Ec.^e^ + - I di((/.„(t)|in-|0„(t)). (1.7) 



Without the driving, Ea = £«, as it should be. By inserting the Fourier expansion (1.5), the 
mean energy takes the form 

{4'a,k\4'a,k) ■ (1-8) 

fe 

This form reveals that the kth Floquet channel yields a contribution + khfl to the mean 
energy, weighted by the Fourier coefficient {<j>a.k\4'a.k)- For the different methods to obtain 
the Floquet states, we refer the reader to the reviews |^,^, and the references therein. 

The invariance of the static Hamiltonian under parity P : {x,p, t) {—x, —p, t) is vi- 
olated by the dipole driving force. With the above choice of the driving, however, a more 
general, dynamical symmetry remains. It is defined by the operation [p" 



PT:{x,p,t)^{-x,-p,t + T/2) (1.9) 

and represents a generalized parity acting in the extended phase space spanned by x, p, and 
phase, i.e., time t mod T. While such a discrete symmetry is of minor importance in classical 
physics, its influence on the quantum mechanical quasispectrum {ea{S, ft)} is profound: It 
devides the Hilbert space in an even and an odd sector, thus allowing for a classification 
of the Floquet states as even or odd. Quasienergies from different symmetry classes may 
intersect, while quasienergies with the same symmetry typically form avoided crossings. The 
fact that Pt acts in the phase space extended by time imodT, results in a particularity: 
If, e.g., \(p{t)) is an even Floquet state, then cxp{int)\ip{t)) is odd, and vice versa. Thus, 
two equivalent Floquet states from neighboring Brillouin zones have opposite generalized 
parity. This means that a classification of the corresponding solutions of the Schrodinger 
equation, \ip{t)) = cxp{—iet / h)\(p{t)) , as even or odd is meaningful only with respect to a 
given Brillouin zone. 



1.3 Coherent tunneling 

With the driving switched off, 5 = 0, the classical phase space generated by T/dw exhibits 
the constituting features of a bistable Hamiltonian system: A separatrix at i? = forms the 
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border between two sets of trajectories: One set, with E < 0, comes in symmetry-related 
pairs, each partner of which oscillates in either one of the two potential minima. The other set 
consists of unpaired, spatially symmetric trajectories, with E > Q, which encircle both wells. 

Torus quantization of the integrable undriven double well implies a simple qualitative pic- 
ture of its eigenstates: The unpaired tori correspond to singlets with positive energy, whereas 
the symmetry-related pairs below the top of the barrier correspond to degenerate pairs of 
eigenstates. Due to the almost harmonic shape of the potential near its minima, neighboring 
pairs are separated in energy approximately by hujQ . Exact quantization, however, predicts 
that the partners of these pairs have small but finite overlap. Therefore, the true eigenstates 
come in doublets, each of which consists of an even and an odd state, |<i>+) and |$~), respec- 
tively. The energies of the nth doublet are separated by a finite tunnel splitting A„. We can 
always choose the global relative phase such that the superpositions 

ieL) = -L(|$+)±|ci,-)) (1.10) 

are localized in the right and the left well, respectively. As time evolves, the states |$^), |$^) 
acquire a relative phase exp(— iA„i/?i) and |<i?5J), |<i?J^) are transformed into one another after 
a time 7r/i/A„. Thus, the particle tunnels forth and back between the wells with a frequency 
A„/fi,. This introduces an additional, purely quantum-mechanical frequency scale, the tun- 
neling rate Ao//i of a particle residing in the ground-state doublet. Typically, tunneling rates 
are extremely small compared to the frequencies of the classical dynamics. 



The driving in the Hamiltonian (1.1), even if its influence on the classical phase space is 
minor, can entail significant consequences for tunneling: It may enlarge the tunnel rate by 
orders of magnitude or even suppress tunneUng altogether. For adiabatically slow driving, i.e. 
$7 <C Ao/?i, tunneling is governed by the instantaneous tunnel splitting, which is always larger 
than its unperturbed value Aq and results in an enhancement of the tunneling rate [^. If the 
driving is faster, the opposite holds true: The relevant time scale is now given by the inverse of 
the quasienergy splitting of the ground-state doublet S/jei — eo|. It has been found that 
in this case, for finite driving amplitudes, |ei— eo| < Aq. Thus tunneling is always decelerated. 
When the quasienergies of the ground-state doublet (which are of different generalized parity) 
intersect as a function of F, the splitting vanishes and tunneling can be brought to a complete 
standstill by the purely coherent influence of the driving — not only stroboscopically, but also 
in continuous time [p|-|rT|]. 

So far, we have considered only driving frequencies much smaller than the frequency scale 
ojo of the relevant classical resonances. In this regime, coherent tunneling is well described 
within a two-state approximation [[TTI]. Near an avoided crossing, level separations may devi- 
ate vastly, in both directions, from the typical tunnel splitting. This is reflected in time-domain 
phenomena ranging from the suppression of tunneling to a strong increase in its rate and to 
complicated quantum beats [O]. Singlet-doublet crossings, in turn, drastically change the 



quasienergy scales and replace the two-level by a three-level structure. 



Three-level crossings 



A doublet which is driven close to resonance with a singlet can be adequately described in 
a three-state Floquet picture. For a quantitative account of such crossings and the associated 
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Figure 1.2: Quasienergies (a) and mean energies (b) found numerically for the driven double well 
potential with D = Es/hujo — 2 and the dimensionless driving strength F = 10^^. Energies of 
states with even (odd) ge neral ized parity are marked by full (broken) lines; bold lines (full and broken) 



correspond to the states (1.16) which are formed from the singlet 
frequency Q > 1.5 ljq corresponds to a detuning S — — E^ - 



) and the doublet 

hQ. < 0. 



A driving 



coherent dynamics, and for later reference in the context of the incoherent dynamics, we shall 
now discuss them in terms of a simple three-state model, which has been discussed in the 
context of chaotic tunneling [^[l^]. In order to illustrate the above three-state model and to 
demonstrate its adequacy, we have num erically studied a singlet-doublet crossing that occurs 
for the double-well potential. Eg. ( |l.l[ ), with D = 2, at a driving frequency w 1.5 t^o and 
an amplitude F = 0.001 (Fig. [L2| ). 

Far outside the crossing, we expect the following situation: There is a doublet (subscript d) 
of Floquet states 

|V.+ (t)>=e-'^^/''|0+(t)), ^^^^^ 

with even (superscript 4-) and odd (— ) generalized parity, respectively, residing on a pair 
of quantizing tori in one of the well regions. We have assumed the quasienergy splitting 
A = e^j^ — (as opposed to the unperturbed splitting) to be positive. The global relative 
phase is chosen such that the superpositions 

|0r.lW> = ^(l0^W)±l0dW» (1-12) 

are localized in the right and the left well, respectively, and tunnel back and forth with a 
frequency A/h. 

As the third player, we introduce a Floquet state 

|V'r(i))=e-(^^+A+W'>|^-(i))^ (1.13) 



6 



1 Decoherence in resonantly driven bistable systems 



located mainly at the top of the barrier (subscript t), so that its time-periodic part 
contains a large number of harmonics. Without loss of generality, its parity is fixed to be 
odd. Note that \(ty^{t)) are in general not eigenstates of the static part of the Hamiltonian, 
but exhibit for sufficiently strong driving already a non-trivial, T-periodic time-dependence. 
For the quasienergy, we assume that — e\ + A. + 5 = e'^ + 5, where the detuning 
5 = EfT — — h£l serves as a measure of the distance from the crossing. The mean energy 
of \^^{t)) lies approximately by hO. above the doublet such that \E'^ — E'^\ ^ E^ — E^. 

In order to model an avoided crossing between |0j ) and \4>^), we suppose that there is a 
non- vanishing fixed matrix element 



1 



T 

di(0dl^Dw|0r) >0. (1.14) 



For the singlet-doublet crossings under study, we typically find that A < 6 ^ hfl. Neg- 
lecting the coupling with all other states, we model the system by the three-state (subscript 3s) 
Floquet Hamiltonian 



Hss = ej + I A b I (1.15) 

b A + 6 

in the three-dimensional Hilbert space spanned by \4'^{t)), \(j)^{t))}. Its Floquet 

states are 

l^lit)) = (|0dW>cos/3-|<^-(O)sin/3), (1.16) 
I02"W> = (|0dW)sin/3+|0rW>cos/3). 
with quasienergies 

and mean energies, neglecting contributions of the matrix element b, 

= ^d' 

= cos^ + sin^ /?, (1.18) 
E:^ = E^ sin^ (3 + E~ cos^ /?. 
The angle (3 describes the mixing between the Floquet states and |(/)^) and is an 



alternative measure of the distance to the avoided crossing. By diagonalizing the matrix (1.15), 
we obtain 

2/3 = arctan ^ , < /? < |. (1.19) 

For /? 7r/2, corresponding to — 5 3> &, we retain the situation far right of the crossing, as 
outlined above, with \4)i) w — 10^), |02") ~ l^d )• To the far left of the crossing, i.e. for 
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/3 ^ or 5 ^ 6, the exact eigenstates and 102^) have interchanged their shape [p[p^. 
Here, we have « and j^j") ~ \(t'^)- The mean energy is essentially determined by 
this shape of the state, so that there is also an exchange of and E2 in an exact crossing, 
cf. Eq. (1.18), while E'^ remains unaffected (Fig. 1.23). 

To study the dynamics of the tunneling process, we focus on the state 



^o+*A|</.+ (i))+c--r*A|<^-(i))cos/3- 



^*/''|02-(i))sin/3). 

(1.20) 



It is constructed such that at t = 0, it corresponds to the decomposition of in the basis 



(1.16) at finite distance from the crossing. Therefore, it is initially localized in the right well 
and follows the time evolution under the Hamiltonian (1.15). From Eqs. (1.12), (1.16), we 
find the probabilities for its evolving into |(/)r), or \4>t), respectively, to be 




COS ^ ^ " ^ COS^ /? + COS — 



4)t 



sin^/? 



cos 



1 — cos 




(1.21) 



cos^ (3 siv? p. 



At sufficient distance from the crossing, there is only little mixing between the doublet and the 
resonant states, i.e., sin /3 ^ 1 or cos /3 <C 1. The tunneling process then follows the familiar 
two-state dynamics involving only \(f>'^) and with tunnel frequency A/h. Close to the 
avoided crossing, cos /3 and sin (3 are of the same order of magnitude, and jf/)^) become 
very similar to one another Each of them has now support at the barrier top and in the well 
region, they are of a hybrid nature. Here, the tunneling involves all the three states and must 
be described at least by a three-level system. The exchange of probability between the two 
well regions proceeds via a "stop-over" at hte top of the barrier 



1.4 Dissipative tunneling 

The small energy scales associated with tunneling make it extremely sensitive to any loss of 
coherence. As a consequence, the symmetry underlying the formation of tunnel doublets is 
generally broken, and an additional energy scale is introduced, the effective finite width at- 
tained by each discrete level. As a consequence, the familiar way tunneling fades away in the 
presence of dissipation on a time scale icoh- In general, this time scale gets shorter for higher 
temperatures, reflecting the growth of the transition rates. However, there exist counterintu- 
itive effects: in the vicinity of an exact crossing of the ground-state doublet, coherence can be 
stabilized with higher temperatures ^ until levels outside the doublet start to play a role. 
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As a measure for the coherence of a quantum system we employ in this work the Renyi 
entropy 

= (1.22) 

In our numerical studies we will use 5*2 which is related to the purity tr(/9^). It possesses a 
convenient physical interpretation: Suppose that p describes an incoherent mixture of n states 
with equal probability, then tr(/5^) reads 1 /n and one accordingly finds 52 = In n. 

Floquet-Markov master equation 



To achieve a microscopic model of dissipation, we couple the driven bistable system ( 1 . 1 ) 
bilinearly to a bath of non-interacting harmonic oscillators ||^,[l^, 16 1. The total Hamiltonian 
of system and bath is then given by 

m - i/ow(.) + 1 + '^^l - ^f)') . 0^23) 

Due to the bilinearity of the system-bath coupling, one can eliminate the bath variables to get 
an exact, closed integro-differential equation for the reduced density matrix p{t) = trBPtotai(i). 
It describes the dynamics of the central system, subject to dissipation. 

In the case of weak coupling, such that the dynamics is predominatly coherent, the re- 
duced density operator obeys in good approximation a Markovian master equation. The Flo- 
quet states {(pait)) form then a well-adapted basis set for a decomposition that allows for an 
efficient numerical treatment. If the spetral density of the bath influence is ohmic [^|l6ll, the 
resulting master equation reads [|l7[ |l8[] 

Paflit) = ~ ^l3)Pa0{t) + ^ Cap ,a' p' Pa' 13' ■ (1.24) 

a'P' 

The time-independent dissipative kernel 

Cap,a'f5' — ^ {Naa',k + Npp'^k) Xaa' ,kXfj' j}^-k 
k 

— Spp' ^ Nfjii^i j.Xal3" ^-kXp"a' ,k (1.25) 

13", k 

— Saa' ^ Na"l3',kXfj'a"-kXa"[3,k 

a"k 

is given by the Fourier coefficients of the position matrix elements, 

Xap.k = Ate-'^'''{Mt)\x\cj,p(t))) = X*p^^_^ (1.26) 

and the coefficients 

Nap,k^N{ta-ep + khn), N{e) = "^7 ee/fceT _ i (1-27) 
which consist basically of the spectral density times the thermal occupation of the bath. 
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t [10'"^ X 2-Kluja] t [lO'"^ X 27r/t^o] 



Figure 1.3: Time evolution of the state at the center of the singlet-doublet crossing found for 
D — 2, F = 10~^, and = 1.5 ljo- The full line depicts the return probability and the broken line the 
occupation probability of the state at the top of the barrier. The dotted line marks the Renyi entropy 52. 
Panel (b) is a blow-up of the marked region on the left of panel (a). 



Dissipative time evolution 



We have studied dissipative tunneling at the particular singlet-doublet crossing introduced in 
Sec. 1.3 (see Fig. pT 3). T he time evolution has been computed numerically by integrating 



the master equation ( |1.24| ). As initial condition, we have chosen the density operator p(0) = 
{4>l\, i-C- a pure state located in the left well. 

In the vicinity of a singlet-doublet crossing, the tunnel splitting increases and during the 
tunneling, the singlet |0t) at the top of the barrier becomes populated periodically with fre- 
quency — £i\/fi, cf. Eq. (1.21) and Fig. 1.3b. The large mean energy of this singlet results 
in an enhanced entropy production at times when it is well populated (dashed and dotted line 
in Fig. L3 b). For the relaxation towards the asymptotic state, also the slower transitions within 
doublets are relevant. Therefore, the corresponding time scale can be much larger than tcoh 
(dotted line in Fig. 1.3a). 

To obtain quantitative estimates for the dissipative time scales, we approximate tcoh by the 
growth of the Renyi entropy, averaged over a time tp, 



1 



t-coh 



1 



1 



P Jo 



dt'S2it') = -{S2{tp)- S2{0}) 



(1.28) 



Because of the stepwise growth of the Renyi entropy (Fig. |l.3| b), we have chosen the propa- 
gation time tp as an n-fold multiple of the duration 27rfi/|e2" — ej^j of a tunnel cycle. For this 
procedure to be meaningful, n should be so large that the Renyi entropy increases substantially 
during the time tp (in our numerical studies from zero to a value of approximately 0.2). We 
find that at the center of the avoided crossing, the decay of coherence, respectively the entropy 
growth, becomes much faster and is essentially independent of temperature (Fig. p~4]a). At a 
temperature fceT = IQ^'^hujQ it is enhanced by three orders of magnitude. This indicates that 
transitions from states with mean energy far above the ground state play a crucial role. 
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Figure 1.4: Decoherence time (a) and Renyi entropy S2 of the asymptotic state (b) in the vicinity of 
the singlet-doublet crossing ior D — 2, F — 10~'^, and f2 = 1.5 ujo- The temperature is given in units 
of hujQ. 



As the dynamics described by the master equation ( 1.24 ) is dissipative, it converges in 
the long-time limit to an asymptotic state Poo{t)- In general, this attractor remains time de- 
pendent but shares the symmetries of the central system, i.e. here, periodicity and generalized 



parity. However, the coefficients ( 1.25| ) of the master equation for the matrix elements are 



time independent and so the asymptotic solution also is. The explicit time dependence of the 
attractor has been effectively eliminated by the use of a Floquet basis. 

To gain some qualitative insight into the asymptotic solution, we focus on the diagonal 
elements 

^aa^a'a' ^ -^aa' ,n l^aa' ,n I ; ^ G; , (1.29) 

n 

of the dissipative kernel. They give the rates of direct transitions from {(pa') to \(f>a)- Within a 
golden rule description, these were the only non-vanishing contributions to the master equa- 
tion to affect the diagonal elements paa of the density matrix. 

In the case of zero driving amplitude, the Floquet states \(pa) reduce to the eigenstates of 
the undriven Hamiltonian Huw- The only non-vanishing Fourier component is then |(/)q,o), 
and the quasienergies reduce to the corresponding eigenenergies Ea- Thus Caa.a'a' only 
consists of a single term proportional to N{Ea — Ea')- It describes two kinds of thermal 
transitions: decay to states with lower energy and, if the energy difference is less than ksT, 
thermal activation to states with higher energy. The ratio of the direct transitions forth and 
back then reads 

^aa,a' a' ( ^a' \ ^r,\ 

— exp ; — — — . (1.30) 



^a'a'^aa \ kj^T 

We have detailed balance and therefore the steady-state solution is 

p„„.(oo)~e-^°/'=^^5„„-. (1.31) 



1.5 Conclusions 
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In particular, the occupation probability decays monotonically with the energy of the eigen- 
states. In the limit k^T 0, the system tends to occupy mainly the ground state. 

For a strong driving, each Floquet state \(j)a) contains a large number of Fourier compo- 
nents and Caa.a'a' IS givcn by a sum over contributions with quasienergies Cq, — ea' + kMl. 
Thus, a decay to states with "higher" quasienergy (recall that quasienergies do not allow for a 
global ordering) becomes possible due to terms with fc < 0. Physically, it amounts to an inco- 
herent transition under absorption of driving-field quanta. Correspondingly, the system tends 
to occupy Floquet states comprising many Fourier components with low index fc. According 
to Eq. (1.8), these states have a low mean energy. 

The effects under study are found for a driving with a frequency of the order ujq. Thus, for 
a quasienergy doublet, not close to a crossing, we have je^ — €a' \ ^ and Ca'a' .aa is dom- 
inated by contributions with n < 0, where the splitting has no significant influence. However, 
except for the tunnel splitting, the two partners in the quasienergy doublet are almost identical. 
Therefore, with respect to dissipation, both should behave similarly. In particular, one expects 
an equal population of the doublets even in the limit of zero temperature in contrast to the 
time-independent case. 

In the vicinity of a singlet-doublet crossing the situation is more subtle. Here, the odd 
partner, say, of the doublet mixes with the singlet, cf . Eq. ( 1 . 1 6 ), and thus acquires components 
with higher energy. Due to the high mean energy of the singlet, close to the top of the 
barrier, the decay back to the ground state can also proceed indirectly via other states with 
mean energy below . Thus, \(f)i) and \ <j)2) are depleted and mainly |0(J") will be populated. 
However, if the temperature is significantly above the splitting 2h at the avoided crossing, 
thermal activation from \<f)^) to \<j)i2)^ accompanied by depletion via the states below , 
becomes possible. Asymptotically, all these states become populated in a cyclic flow. 

In order to characterize the coherence of the asymptotic state, we use again the Renyi en- 
tropy (1.22). According to the above scenario, we expect 5*2 to assume the value In 2, in a 
regime with strong driving but preserved doublet structure, reflecting the incoherent popula- 
tion of the ground-state doublet. In the vicinity of the singlet-doublet crossing where the dou- 
blet structure is dissolved, its value should be of the order unity for temperatures fceT ^ 25 
and much less than unity for k^T ^ 2b (Fig. 1.43). This means that the crossing of the singlet 
with the doublet leads asymptotically to an improvement of coherence if the temperature is be- 
low the splitting of the avoided crossing. For temperatures above the splitting, the coherence 
becomes derogated. This phenomenon compares to chaos-induced coherence or incoherence, 
respectively, found in Ref. for dissipative chaos-assisted tunneling. 



1.5 Conclusions 

For the generic situation of the dissipative quantum dynamics of a particle in a driven double- 
well potential, resonances play a significant role for the loss of coherence. The influence of 
states with higher energy alters the splittings of the doublets and thus the tunneling rates. We 
have studied decoherence in the vicinity of crossings of singlets with tunnel doublets under 
the influence of an environment. As a simple intuitive model to compare against, we have 
constructed a three-state system which in the case of vanishing dissipation, provides a faithful 
description of an isolated singlet-doublet crossing. The center of the crossing is characterized 
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by a strong mixing of the singlet with one state of the tunnel doublet. The high mean energy of 
the singlet introduces additional decay channels to states outside the three-state system. Thus, 
decoherence becomes far more effective and, accordingly, coherent oscillations fade away on 
a much shorter time scale. 
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